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Abstract 

The self-destructive percolation model is defined as follows: Con- 
sider percolation with parameter p > pc- Remove the infinite occupied 
cluster. Finally, give each vertex (or, for bond percolation, each edge) 
that at this stage is vacant, an extra chance 6 to become occupied. Let 
Sc{p) be the minimal value of 5, needed to obtain an infinite occupied 
cluster in the final configuration. This model was introduced some 
years ago by van den Berg and Brouwer. They showed that, for the 
site model on the square lattice (and a few other 2D lattices satisfying 
a special technical condition) that 6c{p) > In particular, 5c{p) 

is at least linear m. p — pc- 

Although the arguments used by van den Berg and Brouwer look 
quite rigid, we show that they can be suitably modified to obtain 
similar linear lower bounds for 5c{p) (with p near pc) for a much larger 
class of 2D lattices, including bond percolation on the square and 
triangular lattices, and site percolation on the star lattice (or matching 
lattice) of the square lattice. 
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1 Introduction 



Some years ago van den Berg and Brouwer, motivated by the study of forest- 
fire processes, introduced the self-destructive percolation model (see [Ij and 
[3]). This model can be described in a few steps as follows. 

Let G = (V, E) be a connected, infinite, locally finite graph. The first 
step is to perform independent site percolation with parameter p on this 
graph: we declare each vertex occupied with probability p and vacant with 
probability 1 —p, independent of the other vertices. U U,W G V, we use the 
notation {U ^ W} for the event that there is a path of occupied vertices 
connecting some vertex of U to some vertex of W, and {U ^ oo} for the 
event that there exists an infinite path of occupied vertices starting from 
some vertex in U. If G = (V, E) is transitive (i.e., without lost of generality 
we can choose any vertex to be the origin), let 9{p) be the probability that 
the origin belongs to an infinite occupied cluster. 

Since 9{p) is monotone in p, there is a critical value Pc G [0, 1] such that 
9{p) > if J9 G {pc, 1] and 9{p) = if p G [0,Pc)- It is well-known (see [1] for 
this and other background results on percolation) that 9{pc) = for a large 
class of graphs including the ci-dimensional hypercubic lattice, for d = 2 or 
d > 19, the triangular and hexagonal lattices. 

Now suppose that, by some 'catastrophe', the infinite occupied clusters 
are destroyed, that is, each vertex that belongs to an infinite occupied cluster 
becomes vacant. After the catastrophe we give each vacant vertex an extra 
chance to become occupied; more precisely, each vacant vertex becomes oc- 
cupied with probability 6, independent of everything else. Let Pp^s be the 
distribution of the final configuration and 9{p, 6) := Pp^s{0 <-> oo). 

An equivalent but more formal (and often more convenient) description 
of the self-destructive percolation model is the following: let Xj, z G V be a 
sequence of i.i.d. 0-1 random variables with parameter p, and let Yi, z G V be 
another sequence of i.i.d. 0-1 random variables with parameter 6. (Here, we 
can interpret as vacant and 1 as occupied). Moreover we take the sequence 
of YiS independent of the sequence of Xj's. Let X*, z G V be defined by 



Here, by X-occupied path we mean a path on which each vertex j has 
Xj = 1. Finally, we define Zi := X* V Yi. This sequence Zi, z G V is the final 




1, if Xj = 1 and there is no infinite X-occupied path from z 
0, otherwise 



(1.1) 



2 



configuration and the measure Pp^ is its distribution. Analogously (with 
obvious modifications of the definitions) one can define the self-destructive 
bond percolation model. 

Monotonicity in 6 obviously implies that for each p G [pc, 1] there exists a 
Sc{p) e [0, 1] such that e{p, 6) = for all 6 G [0, 6cip)) and e{p, 6) > for all 
S G {Sc{p), 1]. It is also clear that if 6 > Pc, then 9{p, 6) > for all p G [0, 1]; 
that is, 

Vp 6,{p) < p,. (1.2) 

As to lower bounds for pc, Proposition 3.1 of [1] and Proposition 2.3.2 in 
[3j show that for site percolation on the square (and the triangular and the 
honeycomb) lattice 5c(p) is at least linear in p — p^. More precisely it was 
shown for these models that 

If p{l -5)> Pc, then e{p, 6) = 0, (1.3) 
which is obviously equivalent to 

Se{p)>^^, P>Pc. (1.4) 
P 

(In fact, it is conjectured in that paper that for self-destructive percolation 
on these lattices Sc{p) does not go to as p | Pc, but a proof or disproof of 
that conjecture seems out of reach at the moment). 

Using totaly different arguments it has been proved that (11.41) also holds 
for site percolation on the binary tree (see Theorem 2.6.1 in [3J or Theorem 5.1 
in pP). This raises the question how general (11.41) . or, at least, the following 
weaker property, are: 

3p>pc3C >0 s.t. Vp G [pcp] Sc{p) > C{p - Pc). (1.5) 

As mentioned in Remark (ii) just after the above mentioned Proposition 
3.1 in [Ij, the argument of this Proposition does not work for bond perco- 
lation on the square lattice. However, in Section 2 we will show that the 
statement of that Proposition is true for that model. We do this by refining 
the argument in [1]. 

On the other hand, it is easy to see that (11.41) can not hold for the bond 
model on the triangular lattice, the site model on the matching (or star) 
lattice of the square lattice, and more generally any percolation model with 
Pc smaller than 1/2: For such models {p — Pc)/p is clearly larger than 1/2 
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(and hence larger than p^) for p sufficiently close to 1. For such p f ll.4l) would 
contradict fll.2l) . 

Nevertheless, still further refinement of the arguments shows that for a 
large class of 2D percolation models, including the above mentioned bond 
model on the triangular lattice and the site model on the star lattice of the 
square lattice, fll.SI) does hold. This is done in Section 3. 

This paper focuses mainly on transitive graphs. However, related to Re- 
mark [6] at the end of Section 3, we make the following comment on non- 
transitive graphs. On such graphs Pp{v oo) may depend on the vertex v, 
and for this reason the notation 6{p; v) is used. However (as is well-known), 
by positive association (FKG) 9{p; w) > Pp{v ^ w)9{p; v) for all vertices 
f , w. Therefore 6{j)] v) is positive if and only if 6{j)] w) is positive, and hence 
Pc does not depend on v. For similar reasons (see for positive association 
for the self-destructive percolation model) 5c{p) does not depend on v either. 

2 Bond percolation on the Square Lattice 

Let = (Z^,E(L^)) be the square lattice, where 

E(L2) = {(x,y);||x-y||i = l} 

and let Gcb = (^^jEc^) be the chess-board lattice, where 

Ec6 = {{x,y); \\x - y\\i = 1} 

U {{ixi,X2), (2/1, ^2)); (2/1,1/2) = (1, 1) + {xi,X2) and Xi + X2 is even } 
U {{{xi,X2), (2/1,2/2)); (2/1,2/2) = (1, -1) + ixi,X2) and Xi + X2 is odd } 

It is well-known that bond percolation on the square lattice is equivalent 
to site percolation on its covering graph (see section 2.5 in [5] or section 1.6 
in [3]), the chess-board lattice. 

The next result is an extension of Proposition 3.1 in [T] for bond percola- 
tion on the square lattice. For the proof of that proposition it was essential 
that the lattice under consideration is a subgraph of its matching lattice. 
(For the notion 'matching lattice', see e.g. Section 3.1 in [1]. In the par- 
ticular case of the square lattice, the matching lattice is obtained from the 
square lattice by adding, in each face, the two diagonals as extra edges; the 
matching lattice of the triangular lattice is the tringular lattice itself: it is 
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self- matching). The chess-board lattice does not have this property. This is 
why the proof of [T] does not work for the site model on that lattice (and 
hence for bond percolation on the square lattice). However, the chess-board 
lattice is a translation ('along an edge') of its matching lattice, and we will 
exploit this property to modify the proof of the above mentioned Proposition 
3.1 in [IJ. 

Theorem 1. For the self-destructive site percolation model on the chess- 
board lattice (or, equivalently , bond percolation on the square lattice), it holds 
that ifp{l-5)> Pc then e{p,5) =0. Hence, 5c{p) > 

Proof. Given v G Z^, we will use the notation v := v -\- (1,0). Recall the 
sequences of 0-1 valued random variables, X^,Y^,Zy v G Z^, introduced in 
section 1. We color each vertex G red if = 1 and Yfj = 0. Then, each 
vertex v will be red with probability p{l — 6), independently of the other 
vertices. Since, p{l — 6) > Pc it follows from ordinary site percolation (see 
section 11.8 in [3]), that, a.s., there is an infinite red cluster, and this cluster 
contains a circuit around the origin. 
Let 7 be such a red circuit. Define 

7:=7+(l,0). 

Note that 7 is a circuit in the matching graph. Let ?) G 7. By construction, 
Yfj = 0. Moreover, since w is a neighbor of the infinite X-occupied cluster, 
X? = 0. Hence Zfj = 0. Summarizing we have that, almost surely, there is 
a Z-vacant circuit in the matching graph which surrounds or contains the 
origin. Hence 9{p,6) = 0. □ 



3 Other 2D lattices 

Note that in the proof of Theorem [1] (as in that of Proposition 3.1 of fT]) the 
definition of red vertices (or edges) was done in such a way that each vertex 
(edge) is red independently of the other vertices (edges). In the current 
Section the color of a vertex will involve the Y values of all its neighbours. 
This strategy, which will be used to show that (11. 5p is true for a large class of 
2D lattices, leads to dependencies which somewhat complicate the analysis. 
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The main result of this section, Theorem [2] below, is stated for three 
well-known lattices, but, as we shall point out in Remark El the result (with 
practically the same proof) holds for a large class of 2D lattices. 

Theorem 2. For self- destructive site percolation on the matching lattice of 
the square or honeycomb lattice, or self- destructive bond percolation on the 
triangular lattice, the following holds: 
There are p > Pc and C > such that 

Mp e [pcP] 5c{p) > C{p - Pc). 

Proof. From now on G = (V, E) denotes the matching lattice of the square 
or honeycomb lattice, or the covering lattice of the triangular lattice, and Pc 
the critical probability for site percolation on G. (See, however. Remark O 
below). 

Recall that for the proofs in Section 2 we introduced a certain colouring 
of the vertices, and that the colours were i.i.d. so that we could compare the 
result of the coulouring with ordinary percolation. In the current situation 
we will again define a colouring, but now the colours are not independent. 
Nevertheless it turns out that we again obtain a suitable comparison with 
ordinary percolation. First some notation. 

For each vertex t> G V, let Dy := {u G V; {v, u) G E}. Let d := \Dy\ (for 
example, ci = 8 for the matching lattice of the square lattice). 

Recall the sequences Xy,Yy, f G V defined in Section 1. Now define the 
sequence of 0-1 random variables Ry, G V as 

^ ^ f 1, if = 1 and = 0, Vm G Dy 
^ \ 0, otherwise 

If Ry = 1 we say that v is i?— occupied (or, simply, that v is red). 

Before we proceed with the proof of the theorem, we first state Observa- 
tion [3] and state and prove Lemma H] below, which will be used later. 

Observation 3. Let j be a circuit in G. Then every path that starts in the 
interior of 7 and ends in the exterior of 7 contains a vertex which has a 
neighbour on 7. 

Lemma 4. Let e > 0. There is a constant such that for the self-destructive 
site percolation model on G with parameters < 6 < pc and p G {pc, 1 — e) 
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the following holds: For every f G V, every finite subset of vertices F C Y 
and every colouring {r^, u & F) of F, 

P{Y, = l;Ru = Tu, ueF) ^ 
P{Y, = 0;Ru = r„, u e F) - ' 

Hence, P(Yy = l\Ru = r^, u E F) < 0^6. 

Proof of Lemma^ If D^, fl F = 0, the random variable is independent 
of the sequence Ru, u E F, and hence 

P(n = 1; Ru = ru,u e F) ^ _J_ 

P{Y, = 0- Ru = r^,ueF) 1-6' 

For the case F' := (1 F ^ ^, we consider two subcases: If there is at 
least one u & F' with Vu = I then obviously 

P(i; = l,i?„ = r„, VmGF) = 0, 

and it is easy to see that P(Yv = 0; Ru = r^, u E F) > 0. 
If r„ = 0, \lu E F' we have that 

P{Y, = l;Ru = Tu, ueF)= P(Y, = l)P{Ru = ru,u e F \Y, = 1) 

< 6P{Ru = ru,ueF\F'\Y, = l) 
= 6P{Ru = ru,ueF\F'), 

and 

P{Y, = 0;Ru = ru,ue F) = P{Y, = 0)P(i?„ = r„,M e F |n = 0) 

> (l-5)P(P„ = r„,MeF; X, = 0,zeD,\Y, = 0) 
= (1 - - pfPiRu = r^,uEF\Y, = 0;X, = 0,zE D,) 
= {l-6){l-pYP{Ru = ru,ueF\F'). 
Combining these two inequahties we have that 

P(n = l;P. = r.,^£P) ^ 5 

P(n = 0; Ru=ru,ueF) " (1 - - p)'^ " 

So the claim of the lemma holds, with the constant q = ti — ^-d■ This 
completes the proof of Lemma^ 
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Now we continue the proof of Theorem [21 Let e > 0. Suppose p G 
{pc, 1 — e), and let 6 be such that p{l — dc^6) > Pc, where is the constant 
given in Lemma |H 

Let V eY, F a finite set of vertices not containing v and r„ G {0, 1}, Vm G 
F. We have 

P{R, = l\R^ = ru,ue F) = P{X, = l;Y, = 0,zeD,\Ru = r„, u G F) 

= P{X, = 1)P{Y, = 0, zeD,\R^ = ru,ueF) 
> p{l — dc^6), 

where in the inequahty we used Lemma|l](and in the equahty the fact that 
is independent of the collection of random variables {Ru, u E F;Yz, z & F>v})- 
As p(l — dc^6) > Pc, the process (i?^, f G V) dominates an i.i.d. process with 
parameter larger than pc- Comparison with ordinary percolation shows that 
a.s. there is an infinite i?-occupied cluster which contains a circuit around 
the origin. Let 7 be such a circuit. Observe that, by the definition of the 
colourings, 7 belongs to an infinite X-occupied cluster. Define 



For each G F we have either X^ = 1, in which case (by the above obser- 
vation) w belongs to an infinite X— open cluster, or we have X^ = 0. In 
both cases X^ = 0. Since also = for all G F, we have = for all 
w G F. By Observation [3] we now conclude that there is no infinite Z— open 
path starting in O. So we have proved that for all e > and all p < 1 — e, it 
holds that 9{p, 5) = if p{l — dc^5) > Pc, that is, 

Sc{p)>^^, P<l-e. (3.7) 

pdCe 

This completes the proof of Theorem O □ 

Remark 5. Note that in fact we have proved something stronger than the 
claim in the theorem, namely that for every e > there is a > such that 
(13.71) holds for all p < 1 — e. 

This can be extended to the result that there is a C > such that Sc{p) > 
C{p — Pc for all p > Pc- Since the in (13.71) goes to 00 as e goes to 0, this 
result does not follow immediately. However, to get the extension, it suffices 
to show that there is an e > such that 6c{p) is bounded away from for 
p > 1 — e. Or, equivalently, that there are p < 1 and 6 > such that 
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9{p, 6) = for all p > p. 



(3.8) 



This can be proved by (e.g.) arguments very similar to those used in the proof 
that (a) implies (i) in Theorem 5.1 in Artem Sapozhnikov (private com- 
munication) has pointed out to us that further refinements of such arguments 
show that, on the d-dimensional cubic lattice, 5c{p) Pc as p —>■ 1. 

Remark 6. TheoremlE holds for a large class of 2D lattices. Essentially we 
only used that for supercritical percolation the infinite cluster a.s. contains 
a circuit around O, which satisfies Observation O This property holds for 
site percolation on the lattices belonging to the family in Theorem 12.1 in 
Kesten's book J^. Informally speaking, this family consists of lattices which 
belong to a pair of matching lattices with certain periodicity and reflection 
symmetry properties (but which are not necessarily transitive). 
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